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Contacting a bilayer ferromagnet with a singlet even-frequency superconductor allows for the real¬ 
ization of an effective triplet odd-frequency superconductor. In this work, we investigate the Joseph¬ 
son effect between superconductors with different symmetries (e.g. odd- versus even-frequency). In 
particular, we study the supercurrent flowing between two triplet odd-frequency superconducting 
leads through a weak singlet even-frequency superconductor. We show that the peculiar temperature 
dependence of the critical current below the superconducting transition of the weak superconductor 
is a signature of the competition between odd/odd-frequency and odd/even-frequency Josephson 
couplings. 

PACS numbers: 74.45.+c, 74.50.+r, 75.70Cn, 74.20.Rp 


I. INTRODUCTION 

It is well known that superconductivity arises from the 
formation of Cooper pairs of electrons, where the wave- 
function of a pair is a function of spin, space and fre¬ 
quency (or time). The Pauli principle tells us that such 
a wavefunction should be antisymmetric. In dirty metals, 
due to multiple scattering events on impurities, the or¬ 
bital part is necessarily symmetric. As a consequence, a 
spin-singlet pairing is even in frequency and a spin-triplet 
pairing is odd in frequency.® In conventional supercon¬ 
ductors (S) the pairing is even in frequency. However, it 
has been predicted that, thanks to the proximity effect, 
one may induce spin-triplet odd-frequency correlations in 
hybrid superconducting/ferromagnetic structures ( S/F ). 
When F is homogeneous, the triplet proximity effect in¬ 
volves electrons of opposite spins and is short-ranged P 
By contrast, an inhomogeneous magnetization, as, e.g., 
in non-collincar bilayer ferromagnets ( F'/F ), also in¬ 
duces long-range triplet correlations between electrons 
with parallel spinsP^ Furthermore, if F' is short whereas 
F is much longer than the coherence length of singlet cor¬ 
relations, “pure” triplet odd-frequency correlations are 
induced at the extremity of the long ferromagnet. Thus 
the S/F'/F structure realizes an effective triplet odd- 
frequency reservoir (St)- hr this paper, we study how to 
probe these odd-frequency correlations. 

Recently long-range supercurrents have been measured 
in trilayer ferromagnetic Josephson junctions, which can 
be viewed as Josephson junctions between two odd- 
frequency reservoirs (St/St junctions).®® While this 
indicates the presence of triplet odd-frequency correla¬ 
tions, the measurements did not present any peculiarities 
as compared to observations made in “classic” Joseph¬ 
son junctions connecting two conventional superconduc¬ 
tors (S/S). Indeed, from a symmetry point of view the 
St/St junction as well as the S/S junction realize a cou¬ 
pling between two reservoirs sharing the same symmetry 


(odd/odd-frequency for St/St and even/even-frequency 
for S/S), yielding similar supercurrent measurements. 
By contrast, the current-phase relation of St/S junc¬ 
tions is predicted to be superharmonicP®3 This speci¬ 
ficity originates from the odd/even-frequency Josephson 
coupling. Namely, the symmetry mismatch between the 
reservoirs prohibits mechanisms involving the transfer of 
a single Cooper pair. Instead, the supercurrent origi¬ 
nates from the coherent flow of an even number of pairs, 
yielding a peculiar 7r-periodic current-phase relation. 

In this work, we explore the competition between 
odd/odd-frequency and odd/even-frequency Josephson 
couplings through the temperature dependence of the 
critical current of hybrid junctions. In particular, we 
study the current through an St/S/St junction, where 
a conventional superconductor of bare critical tempera¬ 
ture T c is sandwiched between two effective triplet odd- 
frequency reservoirs. Such a junction may be realized 
in a S'/F'/F/S/F/F'/S' hybrid junction, where S' are 
conventional superconductors with a critical temperature 
T', see Fig. [I] 

In the following, we assume that T' c T c . Above the 
critical temperature T c of the superconducting layer, an 
effective odd/odd-frequency Josephson coupling builds 
up: the transfer of “odd-frequency” pairs between leads 
happens via virtual Andreev pairs in the island. Be¬ 
low T c , when the central S'-layer is superconducting, 
the quasi-particles above the gap coexist with the even- 
frequency condensate of Cooper pairs. Therefore, an ad¬ 
ditional odd/even-frequency Josephson coupling is gen¬ 
erated at the interfaces between the layer and the leads, 
generating a double St/S Josephson junction. 

We will show that the currents associated with 
odd/odd-frequency and odd/even-frequency Josephson 
couplings are in competition. Besides a peculiar current- 
phase relation, this leads to a suppression of the critical 
current below the transition temperature of the weak su¬ 
perconductor. 

The outline of the paper is as follows: in Sec. [TTJ, we 
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and ipi are the modulus and phase of the superconduct¬ 
ing order parameter in superconducting nodes. The Pauli 
matrices and cr, : (i = x, y , z) act in particle-hole and 
spin space, respectively. Moreover, the spectral current 
between two nodes or leads i and j, 

4 = ^[gj,9i], ( 2 ) 

is related with the normal-state conductance G^ of the 
“connector” between them. The current flowing through 
that connector is 


Iij = -^Im Tr ^4]- 

CcJ>0 


( 3 ) 


Finally, in the case of a superconducting node with 
“bare” critical temperature T ci , the order parameter 
should satisfy the self-consistency equation 


( b) G c„ 

FIG. 1: (a) Realization of an effective triplet odd-frequency 

superconductor St by contacting a ferromagnetic F'/F bi¬ 
layer with noncolinear magnetizations, with a relative angle 
6 = x~x' between their directions, to a singlet even-frequency 
superconductor S', (b) Effective St/S/St double Josephson 
junction formed with triplet odd-frequency superconducting 
leads connected through a singlet even-frequency supercon¬ 
ducting layer. 


ln^ = 2 7 rT^(i-^Tr[r_( 5 + 5 t )]), (4) 

oj> 0 v 1 ' 

where r_ = (r x — ir y )/2. 

Note that Eqs. ([l]) and Q account for the current con¬ 
servation in each node i. namely y~h L t j = 0 is automat¬ 
ically satisfied. 

III. TRIPLET ODD-FREQUENCY RESERVOIRS 


introduce the formalism and, in Sec. Ill we derive the 
Green function of an effective St reservoir. Then, in 
Sec. IV we treat the full St/S/St junction and compute 
both its current-phase relation and its critical current. 
Finally, we briefly discuss metallic junctions in Sec. [Vj 
before concluding in Sec. m 


As mentioned in the Introduction, a trilayer structure 
consisting of a conventional superconductor and two non- 
collinear ferromagnets realizes an effective triplet reser¬ 
voir (St = S'/F'/F), see Fig. 0a). In this section we 
derive the Green function of St- We take the z-axis per¬ 
pendicular to the layers and choose the magnetizations of 
both ferromagnetic layers to lie in the sy-plane, namely 


II. FORMALISM 

Within the quasi-classical theory, the equilibrium 
properties of hybrid superconducting/ferromagnetic 
junctions can be expressed via the quasi-classical Mat- 
subara Green function g , which is a 4 x 4 matrix in the 
particle-hole and spin spaces, and obeys the normaliza¬ 
tion conditions g 2 = 1 and Tr[g] = 0. Within circuit 
theoryp- g takes the value gi in each superconducting or 
ferromagnetic node or reservoir of the circuit. The Green 
functions in the nodes obey the equations 

27rGQ [(a; + hi.cr)T z + A*, 3 *] + ^ 4 = W 

vi 

J 

Here, Si is the mean level spacing in node i , Gq = e 2 /7r is 
the conductance quantum, u> = (2 n + 1)7tT is a positive 
Matsubara frequency at temperature T (n > 0), hi is an 
exchange field acting on the electron spin in ferromag¬ 
netic nodes, and A; = A,(cos <PiT x — siiupiTy), where A,; 


h' = h' (cos x' x + sin x' V ) in F ', ^ 

h = h (cos x x + sin y y) in F. 

Hence £' F = \/D/h' and £p = \/D/h are the ferromag¬ 
netic coherence lengths in F' and F, respectively. 

The first layer F' generates only short-range correla¬ 
tions. Thus, its length l should not exceed the ferromag¬ 
netic coherence length £' F . By contrast, the non-collinear 
second layer F generates triplet correlations with all dif¬ 
ferent spin projections. To filter out only the long-range 
components, its length L needs to be much longer than 

Within the quasi-classical theory, we call gs , 5 _f', and 
gp the Green function in the S' , F', and F layers, re¬ 
spectively. Here, S' is a reservoir. Thus, in the subgap 
regime, oj <C T', the Green function gs takes the form 
gs = cos (pr x — sirup T y = t v , where p is the supercon¬ 
ducting phase. The Green functions in the ferromag¬ 
netic layers will be determined in the following, using 
the formalism introduced in Sec. [TT| In particular, the 
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Green function g F of our effective triplet reservoir is re¬ 
lated to correlations developing at the edge of F, namely 
9t = 9f(L). 

Assuming / < (J,, we can use circuit theory, where the 
F' layer is a ferromagnetic node. Its Green function g F i 
obeys 

[ 2 ^ Q (u + ih'a x ’)T z + * (G s gs + G F g F ),g F -] = 0 , (6) 

where 6 F > is the mean level spacing in F' and Gs (G. F ) 
is the conductance of the S' /F' (F 1 /F) interface. Fur¬ 
thermore, we introduced the short-hand notation a x > = 
cos x' cr x + sin x' &y ■ In the following, we assume that F' 
is more strongly coupled to S', i.e., Gs G Fl and ne¬ 
glect the leakage current at the F'/F interface to obtain 

(uj + ih'a x ’)T z + 7 S^ /-x 

9f' ~ — /, 9 • (') 

+ ih'a x ') 2 + 

Here js = 6 F >Gs/ (2ttGq) is the induced minigap in F'. 
Note that the same Green function with 75 = A > b! = 
Ez describes a superconductor subject to an external 
Zeeman field Ez- The advantage of using an S' /F' bi¬ 
layer is the possibility of realizing both regimes h' < 75 
and b! > 7 s by tuning, e.g., the transparency of the 
S' /F' interface (Gs) or the thickness of the F' layer 
(S F > oc 1 / 0 - 

Having determined g F >, we now turn to the long ferro¬ 
magnetic layer F of length L £ F . Close to the F'/F 
interface, both short- and long-range correlations coex¬ 
ist. The fast oscillatory behavior of the short-range cor¬ 
relations prevents us from directly applying the circuit 
theory. However, within a few £ F from the F'/F in¬ 
terface, the short-range correlations are suppressed, and 
only the non-oscillating long-range triplet correlations 
survive. Then, for £ F -C z -C £/v with = yjD/ 2 fT , 
we find g F (z) ~ g F = const. In particular, if £ F —» 0 , the 
Green function may be considered constant throughout 
the layer. Thus, within a circuit theory approach, the 
long F layer maps to a ferromagnetic node with £ F —> 0 
or, correspondingly, h —> 00. Consequently g F obeys the 
equation, 

[ujt z + iha x r z + 7 F g F ',gr\ = 0 with h —> 00, (8) 

where 7 F = S F G F /( 2 -kGq), and S F is the mean level 
spacing in F. 

To determine g F , we orthogonally decompose the 
Green function g F t with respect to the length scale over 
which it decays in F'. Namely, we write g F ' = g\\ + g_ l, 
where [g^,a x r z ] = 0 and {g±,a x T z } = 0 . Here g\\ con¬ 
tains the long-range correlations, whereas g± contains 
the short-range correlations in F. Then, g\\ may be 
further decomposed by noting that terms oc a x t z can 
be absorbed into h in Eq. Q. Thus, we write g\\ as 

flj| = 9\\ + J°x T z with J = J Tr [ cr x T i5||]- 

In the limit h —> 00, the short-range correlations are 
completely suppressed, while the long-range correlations 


are not affected by h. As shown in appendix [X[ Eq. ([8]) 
may be rewritten in the form 


[ujt z + 7 f <?||, 5 t ] = 0 , ( 9 ) 


where g« (w) = a(u>)r z + iP(oj)a x T lf> can be obtained from 
Eq. Q with 


a(uj) 


/ 3 (w) 


1 (w ± ih') 

2 ^ yj(w±ih ') 2 + 7|’ 



_± 7 s_ 

V( u± ih ') 2 + 7s ’ 


( 10 ) 

( 11 ) 


where 9 = x — x' is the relative angle between the mag¬ 
netization directions of F and F'. 

Finally, the Green function for the effective triplet 
reservoir solving Eq. © reads 

gr = cosh dr. + zsinh i 9 a x r v , (12) 


with 


coshi?(w) 

sinh^(w) 


w + 'y F a(u>) 

y/[u + 7 F a(u })) 2 - 7 f / 3 2 (u)) ’ 
_ 1 fP(u) _ 

\/[ w + 7 F a(w)] 2 - 7 


( 13 ) 

( 14 ) 


The Green function of the effective triplet reservoir is 
thus described by a single angle 1? which depends, how¬ 
ever on all the parameters (h!, 75, xf, 9 ). Note that 
cosh i 9 corresponds to the normal Green function and en¬ 
codes the density of states, whereas sinhz? corresponds 
to the anomalous Green function, describing the induced 
triplet correlations. As /3 oc sind, we see that the triplet 
correlations vanish for collinear F'/F layers (9 = 0 [tt]) 
as expected, while they are maximal for perpendicular 
magnetizations (9 = 7t/2). For simplicity, we will con¬ 
sider only the case 9 = 7r/2 in the following. The gener¬ 
alization to arbitrary angles is straightforward. 

Knowing the Green function of the effective triplet 
reservoir, we can now obtain its density of states (DoS), 

v(e) = i7)Re[coshd(— ie + 0 + )], ( 15 ) 


where is the density of states of the normal metal. As 
the DoS is even in e, we will consider positive energies, 
e > 0, only. 

The functions a(—ie) and / 3 (—ze) possess singularities 
at e = Ef- = h' ± 75 1 , which are inherited by the DoS. 
We will concentrate on the limiting cases h! -C 75 and 
h' 7s, when these singularities are far away from e = 0. 
In particular, for e, 7 p,h' <C 75, we find 


v(e) 






( 16 ) 


Thus, the zero-energy DoS is enhanced as compared to 
the normal state. Furthermore, it displays a broad dip 
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at e = 0. In the opposite regime, for e,7 f,7s •C h', we 
find 




*4) 


1 + 


1 

2 




(17) 


Here as well, the zero-energy DoS is enhanced. However, 
it possesses a narrow peak at e = 0. The enhancement of 
the DoS with respect to its value in the normal state, as 
well as a peak at e = 0 , were discussed in similar models 
of S/F'/F structures with large exchange fields 

Similarly, we may analyze the triplet correlations en¬ 
coded in sinhi?(w). For xf, h! -C 7 s, 


sinh$(w) 




(18) 


Thus, the correlations decay on the energy scale 75 . By 
contrast, for xf <C 7 s <C h', 


sinht?(u;) 


7 s 


h’ 1 -kL 1 I 
n 1 + 7F 1 + 7A 


(19) 


where 7 = 6G/(2itGq). Furthermore, A and (j) are the 
amplitude and the phase of the order parameter in S, 
satisfying Eq. 

We concentrate on the weak coupling regime, 7 <T, 
where it is possible to perform a perturbative expansion 
of g around its bulk value go- To this end, we write 
g = 9o +5i + • ■ • 1 where g\ <C go • Accordingly, the charge 
current may be written in the form Il/r = ^l/r + ^l/r + 
..., where 

4% = ^“ 7rTIm \ Tl [ T z[9i-U 9 l/r]] • (22) 

cd>0 

The bare Green function of the superconducting node S 
reads 


So 


U}T Z + AqT^ 

’ 


(23) 


where Aq(T') solves the standard BCS equation, 



(24) 


In that case, the correlations are reduced as soon asw > 
7 p and then decay more rapidly on the energy scale h!. 


IV. St/S/St JUNCTION 
A. Current-phase relation 


We are now in a position to study the effective 
St/S/St junction presented in the Introduction, see 
Fig. Qb). Within circuit theory, S' is a superconduct¬ 
ing node of bare critical temperature T c and mean level 
spacing S. It is connected to a left and a right effective 
triplet reservoir (St) via connectors of conductances Gl 
and Gr, respectively. Then g, g^, and gR are the Green 
function in the node, the left reservoir, and the right 
reservoir respectively. Here, g l and g r are the Green 
functions derived in the previous section, whereas g will 
be determined in the following. 

For simplicity, here we consider the effective odd- 
frequency triplet reservoirs to be identical by choosing 
i?i(w) = , &r{ui) = $(w) and G = Gl = Gr. However, 
the reservoirs may have different superconducting phases 
fL/R and magnetization axes Xl/r- Without loss of gen¬ 
erality, we choose fL/R = Ff/2 and Xl/r = ±x/2, such 
that f is the phase bias of the junction whereas x is 
the relative angle between the magnetization axes. Then 
9l/r may be written as 


9l/r = coshi?T z -Hsinhi?CT ±x / 2 T ±v ,/ 2 . (20) 
According to Eq. |T]), g obeys 

[ljt, + AT 0 + 7 gL t 9R ’ 9} = ( 21 ) 


while the phase (f> is undetermined for the bare node. 
However, the 17(1) symmetry is broken once the node is 
coupled to the reservoirs. 

Incorporating go in Eq. (22), we obtain A 1 ) = 0. 


Namely, no Josephson coupling exists at first order. In¬ 
deed, due to the symmetry mismatch between the ef¬ 
fective triplet reservoirs and the singlet superconducting 
dot, a single Cooper pair may not carry a current. We 
thus turn to the next order and compute g \. It obeys the 
first order expansion of Eq. (21) supplemented by the 


normalization condition, namely, 


[UJT Z + A O T 0 ,ffi] = -[7 

{go, 91} = 0. 


9L + gR 


AiT 0 , 3 o], (25) 
(26) 


Here Ai is the first order correction to Ao- The system 
is solved by 


1 


9 1 = 




7 

Af L2 


77 (9l+9r) + Ait 0 ,5o 


9o- 


Additionally, the self-consistency equation Q yields 


l ^ "v to cosh, \ / f 

^=-711. ( W 2 +A 2)3/2 J / l 2-, ( W 2 + A 2 )3 /2 




(27) 

Note that Ai/A 0 < 0, i.e., superconductivity is weak¬ 
ened by the coupling. This reduction does not depend on 
the phase bias f, and may be attributed to the gapless 
property of the odd-frequenc y tr iplet reservoirs that was 
discussed in the end of Sec. Ill cf. Eqs. (16) and (FTT 


It features an inverse proximity effect, where the quasi¬ 
particles, existing at zero energy in the leads, weaken 
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superconductivity in the node S. By consequence, the 
effective critical temperature T* of S at finite 7 is de¬ 
creased, 


T; - T ' 

(,) 2 


T r 


(28) 


cj>0 


Note that a dependence of T* on the relative orientation 
of the magnetizations in adjacent layers would arise in 
higher order in 7 . 


Then, incorporating g\ in Eq. (22), we find 


I ^ = — y^ja(T) cosy sin tp 


(29) 


-b(T) [sin tp cos(2(/>) — (cosy + cos p) sin (2^>)] |, 


where 


pc. 


uj> 0 


b(T) = nTj2 


oj> 0 


x/w 2 + A 2 V w 2 + 

sinh 2 i?(w) Aq 
(w 2 + A 2 ) 3 / 2 ' 


(31) 


The first line of (291 may be identified as a quasi-particle 


current between the two triplet reservoirs that does not 
depend on the phase of the central island. By contrast, 
the second line of (29) may be identified as a condensate 


contribution, corresponding to a superharmonic Joseph- 
son effect between the triplet reservoirs and the singlet 
central island, that depends on the phase <f>. 

Current conservation fixes the phase <f> = kn /2 with 
k G Z. Furthermore, energy minimization imposes 


if cos tp + cosy > 0 , 
otherwise. 


(32) 


As a consequence of the odd/even-frequency Josephson 
coupling between two triplet/singlet pairs, <j> is defined 
modulo 7 r instead of 2-7T. 

Inserting Eq. (32) into (29), we obtairP^ 


G 

I^ [a(T) cos x ~ b(T) sign (cos x + cos V 3 )] s i n ( f- 

(33) 

While the quasiparticle contribution is a continuous func¬ 
tion of the phase bias tp, the condensate contribution dis¬ 
plays a jump at cos x + cos P = O. 1 ^ Both a(T), b(T) > 0. 
Thus, when % < w/2 (y > 7 t/2 ), the two contributions are 
opposed for phases tp < n — y (ip < tt — y) whereas they 
have the same direction for phases tp > w — x (p < tt —y). 
Examples of typical current-phase relations are shown in 
Fig. d 


B. Critical current 

To get more insight into the competition between the 
condensate and quasiparticle contributions, we study the 


Idh 



y = 0 


X = n/6 
y = /r/3 
X = n/2 


FIG. 2: Typical current-phase relations. The current is plot in 
units of I c = , yGa(T)/(2e) for dif fere nt magnetization angles 
y and b(T)/a(T) = 0.15, cf. Eq. (33 1 . The critical current is 
achieved at phase tp = 7t/ 2 for x = 0 and 7 t/ 6, and at phase 
tp = 7r — x f°r X = n /3 and n/2. 


critical current I c of the junction, where I c (T,x ) = 
ma x v [I^ 2 >(T,tp,x)]- Based on the considerations above, 
the critical current is achieved for phase bias tp = 7r/2 or 
tp = -ir — X- Namely, I C (T , y) = max[/i,/ 2 ], where 

h(T, y) = |/^(T,7r/2,y)| = |^|a(T)|cosy|-6(T)|, 

h{T, X ) = |/ (2) (T,7r-y,y)| 

7 G 

= — [a(T)| cosy| + 6(T)] siny. 


For a fixed y, as a function of temperature, the criti¬ 
cal current I c lies either on the I\ or on the / 2 branch. 
While the / 2 -branch increases monotonously with de¬ 
creasing temperature, the temperature dependence of the 
Ji-branch is more complicated. Above T*, the /i-branch 
increases as ln(7s/T), for h' ,T <C 7 s, and as (7 f/T) 2 , 
for 7 p -C 7 s,T -C h', with decreasing temperature. At 
T*, it has a cusp. When further decreasing tempera¬ 
ture below T*, it increases much more slowly or even 
decreases. For 7 g h! and T <T*, we find 


/^y) * {| cosy| In -p- 


(34) 


(T* - T) 


(21 cosyj + 1)|, 


which decreases with decreasing temperature for all val¬ 
ues of y. By contrast, for 7 g <C h! and T < T*, we 
find 

h{T ’ x) - £I c,3, (^) 2 (y) 2 { |c “ x| (35) 


-V (r ' r . r) (2|cosxl + l 


)}■ 


where A f = 31£(5)/[7C(3)] 2 ~ 0.5, which slowly increases 
with decreasing temperature for angles y < 7 t/ 3 . 

Which branch the critical current follows is determined 
by the ratio b(T)/a(T). We find that the ratio b(T)/a(T) 
is zero above T* and increases monotonously below T* , 
satisfying b(T)/a(T) < 1. 
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As a consequence, at high temperatures, the critical 
current follows the /i-branch. At lower temperatures, 
one may distinguish two different behaviors depending on 
whether /i(0,x) is larger or smaller than I 2 (0, x)- The 
critical angle \c at which one switches between the two 
cases is given by 

b(0) = |cosxc|(l-sinXc) ^ 

a(0) l + sinxc 


The solution \c of this equation increases from 0 to 7 t/2 
as 6 ( 0 )/a( 0 ) decreases from 1 to 0 , cf. the dependence 
of the r.h.s. of (36 1 as a function of Xc in Fig. [3j For 
angles x < Xc, the critical current lies on the /i-branch 
at all temperatures. By contrast, for angles x > Xc, the 
current switches to the / 2 -branch at the temperature T \2 
determined by 


b(T 12 ) = | cos y |(1 — sin y) 
a(T 12 ) 1 + sinx 


(37) 


Using these considerations, we now consider the tem¬ 
perature dependence of the critical current for the cases 
h' -C 7 s and h' » 75 , assuming 7 ^ -C T* <C maxfys, h'}. 


f(Xc) 



b(T)/a(T ) 



T/T* 


FIG. 4: Plot of b(T)/a(T) as a function of T/T*, for ti, T* <C 
7 s (solid line) and 7 s,T* -C h! (dashed line). 


h! < 7 s and h! 75 , respectively, and for different 
angles X- 


Id I c() 



FIG. 5: Plot of the critical current I c [in units of / c o = 
7G(7FhV7s) 2 /(2e)] as a function of T/T*, for h! , T c * <C 7 s, 
and different angles \ = 0, 7t/6, 7t/4, 7t/3, 27t/5, 7t/2 (from top 
to bottom). 


FIG. 3: Plot of the r.h.s. of (36 1 , /(x=) = |cosXc|(l — 


siiiXc)/(l + sinxc), as a function of Xc- 


V. METALLIC JUNCTION 


For h! <C 7 s, 6(0 )/a (0) decreases as l/ln^s/T*) 
with decreasing T/f'ys- Thus, the critical angle Xc in¬ 
creases. Since I± decreases with decreasing temperature 
below T*, the critical current displays a unique maximum 
at temperatures close to T* for a wide range of angles. 
This non-monotonous temperature dependence provides 
a clear signature of the competition between odd/odd- 
and odd/even-frequency couplings. 

By contrast, for 75 <C h ', we obtain 6(0)/a(0) ~ 1. 
Thus, a finite temperature T\ 2 below which the critical 
current starts rising again rapidly exists for all angles. In 
the intermediate temperature regime T\ 2 < T < T*, the 
current increases very slowly. Though less pronounced 
than in the opposite parameter regime, this peculiar tem¬ 
perature dependence is a signature of the competition 
between the different symmetry couplings. 

The temperature dependence of the critical current de¬ 
scribed above is illustrated in Figs. [5] and [ 6 j for the cases 


While our calculations are limited to tunnel contacts, 
we may generalize our considerations to metallic junc¬ 
tions based on the identification of the different contri¬ 
butions to the current. To do so, we examine the cor¬ 
responding terms in the Josephson energy. The current 
is then obtained by taking a derivative with respect to 
phase. As before, we assume that the angle between 
the magnetizations of adjacent ferromagnets is 9 = 7 r/2, 
cf. H3 

For a metallic S'/F'/F/S/F/F'/S' junction, the 
quasiparticle contribution near the critical temperature 
of the S layer takes the form 

E T = 9 - c qp A r ( J^ e ~ Ls/is cosycos (tp L - tp R ), 

(38) 

where g ~ G/Gq, G is the normal-state conductance of 
the junction, and c qp is a numerical factor of the order 
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Id I c() 

1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 

FIG. 6: Plot of the critical current I c [in units of I c o = 
7 G[7F7s/(fr , A5)] 2 /(2e)], with Aj$ « 1.76T*, as a function 
of T/T*, for 7 f <S 7 s,T* -C h', and different angles x = 
0,7r/6,7r/4, 7 t/3, 27t/5, 7 r/2 (from top to bottom). Note that 
the critical current saturates at T < 7 f (not visible on the 
scale of the figure). 



of unity. Here we write the phases tpL and pn of the left 
and right su percon ductors S' explicitly. Furthermore, L$ 
and is ~ \J D / T * are the length and coherence length of 
the central superconductor, respectively, where D is the 
diffusion coefficient. The first term oc A' describes the 
usual Josephson energy when the central superconduc¬ 
tor is in the normal state. It increases monotonously as 
the temperature decreases below the critical temperature 
T’ c of the leads and saturates at temperatures T <C T'. 
Thus, close to T* -C T', we may neglect its temperature 
dependence.^The second term oc A 2 (T)/T* accounts for 
the reduction the quasiparticle contribution (responsible 
for the triplet supercurrent flow) due to the developing 
of singlet superconducting correlations below T*, when 
A (T) oc y/T* - T6(T* - T ) is finite. 

The condensate contribution consists of three differ¬ 
ent terms. Namely, there is a Josephson coupling be¬ 
tween the left superconductor and the central supercon¬ 
ductor, depending on the phase difference p^ — 0 , as well 
as a Josephson coupling between the central supercon¬ 
ductor and the right superconductor, depending on the 
phase difference 0 — pf j. As the first harmonic in a bi¬ 
layer junction is short-ranged, only the second harmonic 
survives for both of these contributions. Furthermore, 
there is a crossed term, where two pairs from each of the 
outer superconductors recombine in the central super¬ 
conductor. Thus this contribution depends on the phase 
Pl + Pr — 20. It is suppressed with the length of the cen¬ 
tral superconductor on the scale of the coherence length, 
and it depends on the angle between the magnetizations 
of the left and right ferromagnet. As a consequence, the 
condensate contribution takes the form 

E c f nd = -c cond 5 { cos [2(ip L - 0)] (39) 

T cos [ 2 (<j>-ip R )] 

+2e~ Ls ^ s cos x cos(pl + Pr~ 20) j , 


where c cond oc ^f/L^ With tp^ = ~Pr = p/2, the 
expression simplifies to 

£cond = _ 2 c cond ( 7 ^P{cos(p + e _ Wcosx}cos(20). 

(40) 


Then, minimization of the energy with respect to 0 yields 



if cos p + e is As cos x > 0 , 
otherwise. 


(41) 


Finally, the supercurrent accounting for both the 
quasiparticle and condensate contributions reads 


I = -G 


A' - c qp 


A 2 {T) 


e As cos ^ 


2 c cond sign (cos p + e Ls As cos 


sin tp. 
(42) 


Eq. ( |42| has a similar form as the current-phase relation 
Eq. ( |33[ ) in the tunneling regime. For a short S layer 
with L$ <C is-, we find that the critical current is given 
by I(p = 7t/ 2), corresponding to the Ji-branch discussed 
in Sec. |IVB| Thus, the critical current decreases below 
T* as 


I C = G< A'| cosxl - 


A 2 (T) r 


rji*. 


c qp | cosx| + 2 c 1 


cond 


(43) 

(except for angles \ ~ 7 r/ 2 ). As in the tunneling 
regime, this peculiar temperature dependence provides 
a clear signature of the competition between odd/odd- 
and odd/even-frequency couplings. By contrast, in the 
opposite regime, Ls is, the temperature regime be¬ 
low T*, where the Ji-branch dominates, shrinks to zero. 
Thus, the critical current is given by I(p = 7 r /2 + 0 + ), 
corresponding to the / 2 -branch discussed in Sec. |IVB| 
We obtain 


I C = G 


|A'|cosx|e" is/ « s +2c cond 


A 2 (r) ] 

T* /■ 


(44) 


Here, the quasiparticle and condensate contributions add 
up, which leads to an enhancement of I c below T*. 


VI. CONCLUSION 

Using circuit theory, we have proposed a simple model 
for the Green function gx of an effective triplet odd- 
frequency superconducting reservoir St- Then, we have 
studied the coexistence of singlet even-frequency and 
triplet odd-frequency superconducting correlations in an 
St/S/St Josephson junction. 

We predict that the competition between odd/odd- 
frequency and odd/even-frequency Josephson couplings 





















may be observed in a peculiar temperature dependence 
of the critical current of the St/S/St junction below the 
transition temperature T* of the central superconductor. 
For a large range of parameters, the critical current either 
increases very slowly or even decreases when lowering the 
temperature below T*. This is in sharp contrast with a 
conventional S'/S/S' junction, where the superconduct¬ 
ing transition of the central superconductor leads to an 
enhancement of the critical current PS 

We propose to realize such an St/S/St junction by 
fabricating a hybrid S'/F'/F/S/F/F'/S' junction, i.e., 
by inserting a superconducting layer in the middle of 
an S/F'/F/F'/S junction such as the ones presented in 
Refs. EJH 


with g 2 3 4 5 6 7 8 9 10 11 12 = 1 , can be expanded perturbatively in 1/h: 
9 = 9 ^ + (1 /h)g^ + ... Let us now introduce the 
decomposition A = A± + A m, with [A^,a z r z ] = 0 and 

{A±, <t z t z } = 0. Similarly, gW = g jj n) + g^ . 


In the leading order in h, Eq. jAl| ) yields g ( ° ] = 0 , 
while the normalization condition reads (pi 0 ' 1 ) 2 = 1 . 


In the next order, Eq. (Al) yields 
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2tr 2 Tzfl'i 1) + [A, £|j 0) ] = 0, (A2) 


Appendix A: Derivation of the Green function in 
the strong ferromagnet F 

At large h, the solution of the equation 

\ha z r z + A,g\ = 0, (Al) 


while the normalization condition reads {< 7 ^ + 
S|j 1 \ff||° ) } = 0- It is solved with g^ = A||/y^4| (note 

that A| is scalar), g ^ = — (1/2)<t 2 Tz[A_l, and 
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